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ABSTRACT 

A  Kalman  filter  with  assumed  linear  system  dynamics  was 
derived  from  a  spherical  coordinate  state  vector  and  applied 
to  a  fire  control  system  tracking  environment.   Two  models 
of  target  dynamics  were  developed:   a  constant-velocity 
model  and  a  correlated  random  acceleration  model.   The 
predicted  position  performance  of  the  filters  derived  from 
the  two  basic  models  was  optimized  using  Monte  Carlo  methods 
against  a  set  of  test  trajectories  that  include  constant- 
velocity  maneuvering  target  profiles.   The  results  of  the 
Monte  Carlo  simulation  of  these  filters  are  compared  with  the 
results  obtained  from  a  filter  derived  from  a  Cartesian 
Coordinate  state  vector.   Switch-on-range  adaptive  filters 
were  developed  from  the  two  basic  models  and  evaluated  by 
Monte  Carlo  methods.   The  results  of  the  simulations  of  the 
adaptive  filters  with  the  two  different  state  vectors  are 
compared . 
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I.   INTRODUCTION 

A.   AN  ALTERNATIVE  TO  NONLINEAR  ESTIMATION  IN  GUN  FIRE 
CONTROL  SYSTEMS 

The  increased  performance  characteristics  of  modern 
aircraft  and  missiles  has  made  apparent  the  need  for  more 
effective  fire  control  systems.   With  the  introduction  of 
digital  fire  control  computers,  Kalman  filtering  tech- 
niques offer  an  attractive  means  of  improving  the  effective- 
ness of  a  fire  control  system.   Indeed,  several  recent 
studies  conducted  at  the  Naval  Postgraduate  School  [Refs. 
1-6]  have  specifically  addressed  the  application  of  Kalman 
filtering  techniques  to  fire  control  systems.   It  is 
interesting  to  note  that  the  references  cited  above  all  have 
one  distinctive  feature  in  common;  they  use  Cartesian 
Coordinate  state  vectors.   A  Cartesian  Coordinate  state 
vector  offers  the  attraction  of  linear  plant  dynamics  for  a 
constant-velocity  target  but  has  the  disadvantage  of  non- 
linear measurement  equations  when  used  in  a  fire  control 
system  that  measures  the  target's  position  in  spherical 
coordinates.   One  possible  solution  to  the  nonlinear  measure 
ment  problem  is  to  process  the  sensor  measurements  through 
a  prefilter  as  in  Figure  1.   The  output  of  the  sensor 
consists  of  noise  corrupted  measurements  of  the  target's 
position  in  spherical  coordinates 
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s1(K)  =  r(K)  +  n]L(K) 

s2(K)  ■  (J>(K)  +  n2(K)  (1.1) 

s3(K)  =  6(K)  +  n3(K) 

which  serve  as  inputs  to  the  prefilter.   The  appropriate 
coordinate  transformation  equations  from  Figure  2  are 

x(K)  =  r(K)sin({)(K)cose(K) 

y(K)  =  r(K)coscJ)(K)cose(K)  (1.2) 

z(K)  =  r(K)sin6(K) 

where  it  has  been  shown  [Ref.  3]  that  v(K),  the  measurement 
noise  process  presented  to  the  tracking  filter  is  a  non- 
linear, position-dependent  function  of  the  variances  of 
the  range,  bearing,  and  elevation  errors.   Clearly  then,  the 
measurements  become  correlated  in  the  coordinate  conversion 
process.   Full  implementation  of  such  a  filtering  scheme 
could  conceivably  involve  much  on-line  computation  since 
the  statistics  for  v(K)  are  no  longer  stationary.   Another 
distinct  disadvantage  of  the  prefiltering  method  lies  in 
the  errors  generated  due  to  the  added  finite  precision 
arithmetic  of  the  coordinate  conversion  process.   Such 
errors  may  be  qualitatively  described  as  a  noise  process 
whose  statistical  properties  are  difficult  at  best  to 
identify.   An  excellent  treatment  of  this  specific  topic  is 
found  in  Gelb  [Ref.  7] . 
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Another  approach  to  the  nonlinear  measurement  process 
is  to  linearize  the  measurement  equations  about  some 
reference  or  nominal  trajectory.   Unfortunately,  such  an 
approach  is  usually  accomplished  by  the  implementation  of 
an  extended  or  iterative  extended  Kalman  filter  with  the 
attendant  increase  in  on-line  computing  requirements. 

Consider  the  filtering  problem  with  a  spherical 
coordinate  state  vector 


x(K)  =  [r,r,r,(t,,(i),i,0,0,e]T  (1.4) 


Such  a  state  vector  has  the  advantage  of  linear  measurement 
equations  eliminating  the  correlation  and  position  dependence 
of  the  measurement  noise  process.   The  disadvantage  of  a 
spherical  state  vector  lies  in  the  coupled  and  nonlinear 
system- target-dynamics ,  even  for  a  constant-velocity  target. 
Once  again,  the  implementation  of  this  filter  would  normally 
be  with  an  extended  or  iterative  extended  Kalman  filter. 
This  thesis  proposes  a  filter  having  a  spherical  coordinate 
state  vector  with  assumed  linear  system  dynamics.   In  view 
of  the  fact  that  any  system  model  is  an  approximation  to  a 
physical  system,  this  seems  to  be  a  reasonable  assumption. 
Implementation  of  this  would  be  a  great  deal  simpler  than  the 
previously  mentioned  filters.   A  filter  designed  on  this 
basis  would  then  of  necessity  be  sub-optimal.   However, 
given  acceptable  performance,  the  sub-optimal  filter  might 
be  preferred  to  the  optimal  nonlinear  filter.   Faced  then, 
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with  a  sub-optimal  filter  due  to  modeling  inaccuracies, 
a  technique  must  be  found  to  compensate  for  the  unmodeled 
nonlinear  properties  of  the  system.   Jazwinski  [Ref.  8]  has 
proposed  two  methods  for  the  compensation  of  inaccurate 
a  priori  knowledge  of  the  system.   The  first  involves  a 
fictitious  noise  input  to  the  system,  while  the  second 
overweights  the  most  recent  data.   The  net  result  of  either 
method  is  an  increase  in  the  filter  bandwidth,  making  the 
filter  more  responsive  to  system  transients  but  increasing 
its  sensitivity  to  the  measurement  noise  process.   The 
fictitious  noise  input  method  was  chosen  as  the  most  promis- 
ing and  was  used  in  this  thesis.   So  that  some  quantifiable 
measure  of  the  feasibility  of  these  filters  is  possible, 
this  thesis  paralleled  Parr's  work  [Ref.  4]  and  comparisons 
are  made  with  his  results. 

B.   LINEAR  FILTERING  THEORY 

The  system  studied  in  this  thesis  is  composed  of  a 
target  and  a  sensor  which  gives  noise-corrupted  measure- 
ments of  the  target's  position  in  spherical  coordinates. 
The  system  may  be  approximated  by  the  linear  difference 
equations 

x(K+l)  =  <|>(K+l,K)x(K)  +  r(K)w(K) 

(1.5) 
z(K)  =  H(K)x(K)  +  v(K) 

The  statistics  assumed  for  the  stochastic  processes  w(K) 
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and  v(K)  are 

E[w(K)]  =  E[y(K)]  =  0    for  all  K  (1.6) 

E[v(k)vT(j)  =  R(K)6kj  (1.7) 

E[w(k)wT(j)  =  Q(K)6kj  (1.8) 

E[w(K)vT(K)]  =  E[v(K)wT(K)]  =  0   for  all  K  (1.9) 

E[w(K)xT(0)]  =  E[x(0)wT(K)]  =  0   for  all  K  (1.10) 

E[v(K)xT(0)]  =  E[x(0)vT(K)]  =  0   for  all  K  (1.11) 

The  linear,  discrete  Kalman  filter  is  characterized  by  the 
difference  equations 

x(K/K)  =  x(K/K-l)  +  G(K)[z(K)  -  H(K)x(K/K- 1) ]  (1.12) 

with  the  one-step  prediction  given  by 

x(K/K-l)  =  <KK,K-l)x(K-l/K-l)  (1.13) 

The  Kalman  Gain  matrix,  G(K),  is  found  by  solving  the 
following  equations 

G(K)  =  P(K/K-1)HT(K) [H(K)P(K/K-1)HT(K)  +  R(K)]"1  (1.14) 

P(K/K)  =  [I  -  G(K)H(K)]P(K/K-1)  (1.15) 


P(K+1/K)  =  4)(K+1,K)P(K/K)4>T(K+1,K)  + 


r(K)Q(K)rT(K) 


(1.16) 


To  start  the  gain  calculations,  the  matrix  P(K/K-1)  must 
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be  initialized.   Taking  the  standard  approach,  P(0/-1) 
is  given  by 

P(0/-1)  =  E[(x(0)  -  xQ)(x(0)  -  x0)T]         (1.17) 

where 

E[x(0)]  =  xQ  (1.18) 

An  unbiased  estimator  is  obtained  if  the  filter  is  started 
by  setting 


x(0/-l)  =  xQ  (1.19) 
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II.   DEVELOPMENT  OF  THE  FILTER  MODELS 

A.   CONSTANT-VELOCITY  (1/S2)  MODEL 

The  first  filter  tested  was  based  on  a  dynamic  system 
model  that  assumes  linear  constant-velocity  target  motion 
in  each  of  three  coordinate  directions. 

The  state  transition  matrix  <f>(K+l,K)  and  the  random 
forcing  input  distribution  matrix  r  (K)  are 


<KK+1,K) 


1 

T 

0 

0 

0 

0 

0 
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T 
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1 
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0 
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0 

0 

1 

T 

0 

0 

0 

0 

0 

1 

(2.1) 


and 


r(K) 


0 

0 

0 

0 

rp  Lt 

2 

0 

T 

0 

0 

T2 
2 

0 

T 

(2.2) 


where  T  is  the  sampling  period. 
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Since  the  state  vector  is  expressed  in  spherical 
coordinates,  the  measurement  process  is  linear  and  is  given 


by 


z(K)  =  H(K)x(K)  +  v(K) 


(2.3) 


The  measurement  noise  processes  for  the  individual 
coordinate  directions  are  assumed  to  be  independent 
processes.   The  covariance  of  measurement  noise  matrix 
R(K)  is  diagonal  and  is  given  by 


R(K)  = 


R 


a 


B 


(2.4) 


Because  of  its  simplicity  the  constant-velocity  model 
cannot  realistically  be  expected  to  perform  well  against 
all  types  of  target  trajectories  and  may  not  perform  accep- 
tably against  maneuvering  targets.   The  performance  of  the 
constant-velocity  model  filter  can  be  improved  by  the  use  of 
the  system  random  input  covariance  matrix  Q(K).   Indeed, 
one  of  the  functions  of  the  Q(K)  matrix  is  to  compensate 
for  the  inadequacies  of  a  simplified  model,  and  for  this 
specific  thesis,  Q(K)  can  also  partially  account  for  the 
system  nonlinearities  not  explicitly  included  in  the 
system  model.   The  subject  of  the  sensitivity  of  the 
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filter's  N-step  prediction  performance  to  the  Q(K)  matrix 
forms  the  basis  for  the  simulations  of  Chapter  III. 

Initialization  of  the  filter  requires  a  priori  knowledge 
of  the  statistics  of  the  system  initial  conditions  so  that 
the  gain  calculations  may  be  started.   For  any  fixed  coordi- 
nate system,  it  is  reasonable  to  consider  the  mean  of  all 
possible  system  initial  conditions  as  zero  or 

E[x(0)]  =  xQ  =  0  (2.5) 

It  has  been  shown  [Refs.  7  and  9]  that  under  proper 
initialization,  the  filter  will  remain  unbiased  for 
K  >  0,  i.e. , 

E[e(K/K)]  =  0    for  all  K  (2.6) 

where 

:-  e(K/K)  =  x(K/K)  -  x(K)  (2.7) 

The  unbiased  condition  is  satisfied  if 

x(0/-l)  -  xQ  (2.8) 

Initialization  of  the  covariance  of  estimation  error 
matrix  requires  a  somewhat  different  approach.   Following 
Parr  [Ref.  4]  it  is  quantitatively  shown  that  if  the 
diagonal  elements  of  the  initial  covariance  of  estimation 
error  matrix  are  made  very  large,  the  net  effect  is  that 
the  first  position  estimate  is  made  identically  equal  to 
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the  first  system  measurement.   In  this  manner  the  system 
is  forced  into  tracking  the  target  rapidly.   For  all  simu- 
lations conducted  then,  it  was  assumed  that 


P(0/-1)  = 


(2.9) 


B.   CORRELATED  RANDOM  ACCELERATION  MODEL 

The  constant-velocity  model  discussed  in  the  previous 
section  views  any  target  maneuvers  as  small  perturbations 
about  a  nominal,  constant-velocity  trajectory.   This  tacitly 
assumes  that  any  target  maneuvers  are  uncorrelated  in  time, 
or  from  one  measurement  to  the  next.   Target  maneuvers  will 
generally  take  place  over  a  period  of  n  measurements  and  the 
magnitude  of  n  is  dependent  on  the  level  of  the  target's 
maneuver,  i.e.,  for  a  slow  turn  n  is  relatively  large 
whereas  for  an  evasive  type  maneuver  n  is  smaller.   The 
fact  that  the  accelerations  occur  over  a  finite  period  of 
time  correlates  the  target's  accelerations  (maneuvers)  in 
time.   One  approach  to  modeling  the  effects  of  time- 
correlated  target  maneuvers  is  to  use  a  constant-velocity 
plant  driven  by  correlated  noise  obtained  from  a  low-pass 
filter  of  the  form 


T^  =  IF^7 


(2.10) 
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driven  by  white  noise.   The  constant  a  is  interpreted  as 
the  time  constant  or  degree  of  time  correlation  of  the 
target's  maneuver. 

For  one  dimension,  the  overall  model  transfer  function 
is 


T(s)  = 


a 


s  (s  +  a) 


(2.11) 


as  depicted  in  Figure  3. 

For  the  three-dimensional  case,  the  filter  is  ninth 
order  and  has  the  state  vector 


•   *•       •   ••       *   • •   rp 

x(K)  =  [r,r,r,<f>,<J>,<|>,e,e,e] 


(2.12) 


The  correlated  random  acceleration  model  is  discussed  in 
detail  by  McKinley  [Ref .  5]  and  Singer  [Ref .  10] .   From  the 
cited  references,  the  (3  x  3)  state  transition  matrix 

2.(T)  and  the  (3  x  1)  random  forcing  input  vector  T . (T) 
are  for  the  single  dimension  with  assumed  linear  system 
dynamics 


4>i(T) 
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and 


r.(T) 


T*   T   1   rn     -aT> 

2   "  a  +  —   (1  "  e   > 

a 


T  -  I  (1  -  e'aT) 


(2.14) 


(1  -  e   ) 


For  the  ninth-order  filter  the  state  transition  matrix 
<f>(K+l,K)  is  formed  by 


*(K+1,K)  = 


Jr(T)     ? 


MT)    o 


:b 


$E(T) 


(2.15) 


where  <f>  (T)  ,  <J>g(T),  and  <J>E(T)  are  individually  given  by 
({>•  (T)  in  Equation  (2.13).   In  a  similar  manner,  for  the 
ninth-order  filter  the  three-dimensional  system  matrix 
T (K)  is  given  by 


r(K)  = 


r  (t) 

~r  v    J 

0 

0 

0 

Eb(t> 

0 

0 

0 

rE(T) 

(2.16) 


where  r  (T) ,  r  (T),'T  (T)  are  individually  given  by  r. (T) 
in  Equation  (2.14). 
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The  covariance  of  measurement  noise  matrix  R(K)  is  the 
same  as  given  in  the  preceding  section  by  Equation  (2.4). 
Initialization  of  the  filter  is  carried  out  in  exactly  the 
same  fashion  as  described  in  the  previous  section  by 
Equations  (2.5),  (2.8),  and  (2.9).   As  in  the  preceding 
section,  Q(K),  the  covariance  of  the  random  forcing  input 
white  noise  process  presented  to  the  coloring  filter,  was 
used  to  compensate  for  the  unmodeled  dynamic  nonlinearities . 
The  sensitivity  of  the  filter's  N-step  prediction  perfor- 
mance to  the  Q(K)  matrix  selected  is  treated  in 
Chapter  III. 

In  the  interest  of  saving  computer  time,  the  correlated 
random  acceleration  model  was  tested  using  a  single  value 
for  a,  the  maneuver  time  constant.   Since  Parr  [Ref.  4] 
reported  that  the  best  overall  performance  in  his  simula- 
tions was  obtained  from  a  value  of  a  =  0.5,  that  value  was 
used  in  the  study  of  this  model.   Additionally,  the  maneuver 
time  constant  a,  was  fixed  at  the  same  value  for  each  of 
the  three  coordinate  directions,  though  this  constraint 
need  not  necessarily  be  imposed  on  the  system. 
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III.   RESULTS  OF  THE  MONTE  CARLO  SIMULATION 
OF  THE  PROPOSED  LINEAR  FILTERS 

A.   EXPERIMENTAL  PROCEDURE 

Monte  Carlo  simulation  techniques  were  used  in  the 
optimization  and  testing  of  the  filters  developed  in 
Chapter  II.   A  listing  and  description  of  the  FORTRAN 
program  used  is  given  in  Appendix  A.   All  results  are  the 
products  of  100-member  ensemble  simulations.   Test  trajec- 
tories 1-8  were  used  in  the  optimization  and  testing  process 
and  Appendix  B  gives  a  description  of  the  profile  for  each 
trajectory.   The  profiles  of  trajectories  1,  2,  3,  and  4  are 
identical  to  the  short-range  constant-velocity,  short-range 
maneuvering,  long-range  constant-velocity,  and  long-range 
maneuvering  profiles,  respectively,  used  by  Parr  [Ref.  4]. 
Trajectories  5-8  were  included  to  examine  the  filter's 
performance  against  more  severe  target  maneuvers. 

The  elements  of  R(K) ,  the  covariance  of  measurement 
noise  process  matrix,  were  fixed  for  all  simulations  at 

aR  =  3.0  yds . 

aR  =  0.01722  Deg.  (3.1) 

aE  =  0.01722  Deg. 

The  values  given  in  Equation  (3.1)  are  those  used  by  Parr 
[Ref.  4]  and  are  used  in  this  thesis  in  an  attempt  to 
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duplicate  Parr's  testing  environment,  and  minimize  procedural 
discrepancies . 

The  performance  evaluation  of  the  filters  developed  in 
Chapter  II  is  based  on  their  N-step  prediction  or  positicr.- 
state  prediction  accuracy.   A  detailed  description  of  the 
method  used  in  the  determination  of  the  error  in  the 
predicted  position  is  given  in  Appendix  C.   The  results 
presented  in  this  chapter  are  the  percentages  of  "hits11 
and  "scares"  computed  over  the  entire  ensemble.   A  "hit" 
is  defined  as  an  error  in  the  predicted  position  of  less 
than  eight  yards,  while  a  "scare"  is  defined  as  an  error 
in  the  predicted  position  of  less  than  16  yards,  but  greater 
than  eight  yards . 

Initialization  of  each  filter  simulated  "-as  identical 
in  all  trials  and  is  described  in  Chapter  II  by  Equaticr. ; 
(2.5)  ,  (2.8)  ,  and  (2.9) . 

The  elements  of  the  covariance  of  the  random  forcing 
input  matrix  Q(K)  were  determined  by  a  tv;o- dimensional 
optimization  search.   The  assumption  of  linear  uncoupled 
dynamics  eliminates  any  spatial  correlation,  resulting  in 
a  diagonal  Q(K)  matrix  of  the  fom 


Q(K)  = 


~qr 

0 

0 

0 

% 

0 

0 

0 

q 

(3.2) 
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If  the  system  dynamics  were  truly  linear,  then  the  diagonal 
elements  of  (3.2)  would  represent  the  variances  of  the 
random  forcing  input  to  the  system.   However,  since  Q(K) 
is  also  used  to  compensate  for  inaccuracies  in  the  system 
model,  q  ,  q™  and  qp  must  be  adjusted  accordingly.   To 
simplify  the  optimization  process  for  both  models,  qR  and  qp 
were  constrained  to  be  equal  throughout.   The  two-dimensional 
search  optimized  the  filter's  position-state  prediction 
accuracy  as  a  function  of  two  variables,  q  and  the  pair 

qB  ■  qE. 

In  the  case  of  the  constant-velocity  model,  the  elements 

of  the  matrix  Q(K)  serve  as  a  means  by  which  the  bandwidth 

of  the  filter  is  increased  to  ensure  that  the  filter  will 

respond  to  target  accelerations.   This,  of  course,  is  in 

addition  to  compensation  for  the  unmodeled  nonlinearities  of 

the  system  dynamics.   Because  of  the  dual  role  played  by  the 

matrix  Q(K),  a  brute  force  search  was  done  to  determine  the 

elements  of  Q(K)  that  gave  the  best  position-state  prediction 

performance.   In  a  similar  manner,  for  the  filter  based  on 

the  correlated  random  acceleration  model,  the  bandwidth 

must  be  increased  to  permit  sufficient  response  to  targets 

exhibiting  non-zero  acceleration  rates.   Here  again  brute 

force  search  techniques  were  used  to  determine  the  Q(K) 

matrix  that  gave  the  best  position  prediction  performance. 

For  both  models,  once  the  optimum  Q(K)  matrix  was 

identified  for  a  specific  trajectory,  that  particular 
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filter  configuration  was  tested  against  the  remaining 
trajectories.   The  net  result  being  that  each  optimized 
filter  was  tested  against  seven  additional  trajectories  to 
ensure  that  the  optimum  Q(K)  had  been  found,  and  to  deter- 
mine the  effect  of  using  a  Q(K)  selected  for  a  track  of 
one  type  against  targets  having  different  maneuvering 
characteristics . 

B.   SIMULATION  RESULTS  FOR  THE  CONSTANT-VELOCITY  (L/S2)  MODEL 

The  results  presented  in  this  section  are  for  the  opti- 
mized filter  configurations.   The  term  "Filter  1"  denotes 
the  filter  with  the  specific  Q(K)  matrix  that  gave  the  best 
position  prediction  performance   against  trajectory  1.   Thus, 
the  simulation  results  for  eight  filters,  or  more  precisely, 
eight  filter  configurations  are  presented. 

As  pointed  out  in  the  preceding  section,  the  assumption 
of  linear  uncoupled  dynamics  results  in  the  diagonal  Q(K) 
matrix  given  by  Equation  (3.2).   Table  I  gives  the  diagonal 
elements  of  the  Q(K)  matrix  that  provided  the  best  position 
prediction  performance  for  each  test  trajectory.   The  rela- 
tive magnitude  of  the  diagonal  elements  for  the  optimal 
filters  are  as  expected.   That  is  to  say,  the  elements  of  Q(K) 
for  the  non-maneuvering  trajectories  (1,  3,  7,  8)  are  lower 
in  magnitude  than  those  for  the  maneuvering  trajectories 
(2,  4,  5,  6).   Similarly,  the  elements  of  Q(K)  for  the  long- 
range  maneuvering  trajectories  are  lower  in  magnitude  than 
those  for  the  short-range  maneuvering  trajectories.   Similar 
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Table  I 
DIAGONAL  ELEMENTS  OF  THE  OPTIMAL  Q(K)  MATRICES 
FOR  THE  CONSTANT-VELOCITY  (l/S2)  MODEL 


FILTER  qn  q22  q33 

1        1.00  x  10"6      8.58  x  10"9      8.58  x  10"9 


0.05  0.1153  0.1153 


1.00  x  10'6      1.13  x  10'10     1.13  x  10"10 


0.025  0.006553         0.006553 


0.0075  0.1372  0.1372 


1.00  x  10"4      5.55  x  10"7      5.55  x  10_7 


0.01  3.37  x  10"7      3.37  x  10_7 


0.005  1.66  x  10~7      1.66  x  10_7 
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observations  may  be  made  for  the  moderate  level  maneuvers 
as  opposed  to  the  evasive  level  maneuvers.   Note  that  the 
elements  of  Q(K)  for  the  long-range  evasive  maneuver  profile 
(trajectory  6)  are  smaller  than  those  for  the  long-range 
moderate  level  maneuver  profile  (trajectory  4).   This  appears 
to  be  a  reasonable  phenomenon  in  light  of  the  fact  that 
overweighting  the  estimates  for  a  faster  more  highly  maneu- 
vering target  at  long  range  could  result  in  larger  errors 
in  the  predicted  position.   It  follows  then  that  the  band- 
width of  the  optimal  filter  for  the  long-range  evasive 
maneuver  should  be  somewhat  smaller  than  the  bandwidth  for 
the  moderate  level  maneuver  at  long  range. 

Table  II  gives  the  simulation  results  for  the  optimal 
filter  configurations  as  well  as  for  each  optimal  filter 
configuration  against  the  other  trajectories.   The  results 
shown  are  not  surprising  in  that  for  a  specific  trajectory, 
the  optimal  filter  configuration  gave  the  best  results.   In 
a  few  isolated  cases,  other  filter  configurations  gave 
results  that  approach  those  for  the  optimal  configuration. 
Though  this  was  not  expected,  the  differences  are  not  con- 
sidered small  enough  to  justify  any  additional  effort  in 
the  optimization  process. 

Table  III  compares  the  hit/scare  percentages  from  the 
simulations  of  the  two  models  using  the  two  different  state 
vectors.   The  results  for  the  spherical  coordinate  state 
vector  represent  the  optimal  performance  against  a  specific 
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trajectory.   The  results  tabulated  for  the  Cartesian 
coordinate  state  vector  are  the  best  of  the  results  reported 
by  Parr  [Ref .  4] . 

Tables  IV  and  V  show  the  tracking  performance  of  each 
optimal  filter  against  its  specific  trajectory.   For 
example,  the  tracking  performance  of  filter  1  against  tra- 
jectory 1  is  given,  and  so  on.   The  tabulated  data  represent 
the  means  of  estimation  error,  defined  by  Equation  (2.7),  time 
averaged  over  the  entire  trajectory,  and  the  RMS  value  of  the 
deviation  of  the  means  of  estimation  error  about  their  time 
average. 

C.   SIMULATION  RESULTS  FOR  THE  CORRELATED  RANDOM 
ACCELERATION  MODEL 

The  results  presented  in  this  section  are  for  the  optimized 
filter  configurations.   Since  the  optimization  search  for  the 
ninth-order  filter  required  a  great  deal  more  computer  time, 
this  model  was  only  optimized  against  trajectories  1-3.   As 
in  the  preceding  section,  the  term  "Filter  1A"  denotes  the 
filter  with  the  specific  Q(K)  matrix  that  gave  the  best 
position  prediction  performance  against  trajectory  1. 

The  tabulated  data  of  Table  VI  gives  the  elements  of 
the  Q(K)  matrices  that  provided  the  optimum  performance 
against  each  test  trajectory.   The  data  appear  to  be  con- 
sistent within  the  tracking  environments  presented  to  the 
filter.   Notice  that  in  each  configuration  the  elements  of 
the  optimum  Q(K)  matrices  for  the  correlated  random 
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Table  VI 

DIAGONAL  ELEMENTS  OF  THE  OPTIMIZED 
Q(K)  MATRICES  FOR  THE  CORRELATED 
RANDOM  ACCELERATION  MODEL  FILTERS 


FILTER 

qll 

q22 

q33 

1A 

l.OOxlO"7 

5.00xl0~9 

5.00x10 

2A 

0.100 

14.599 

14.599 

3A 

l.OOxlO"3 

1.46xl0~7 

1.46x10 
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acceleration  model  are  larger  in  magnitude  than  those  for 
the  corresponding  entries  of  the  constant-velocity  model. 
This  seems  reasonable  because  this  filter  is  required  to 
respond  to  a  transient  (due  to  the  non-zero  acceleration 
rate)  that  is  expected  to  be  smaller  in  magnitude  than  the 
transient  encountered  by  the  constant-velocity  model  filter. 
For  this  reason  the  bandwidth  of  this  filter  must  be  in- 
creased to  ensure  sufficient  response  to  the  non-zero 
acceleration  rates.   During  the  optimization  search,  it  was 
observed  that,  with  one  exception,  this  filter  model  is 
extremely  sensitive  to  variations  of  elements  q?~   and  q__ 
corresponding  to  acceleration  rates  normal  to  the  Line-of- 
Sight.   This  is  to  be  expected,  since  such  accelerations 
cause  the  greatest  displacement  from  the  nominal  trajectory 
adding  a  source  for  potentially  large  errors  in  the  pre- 
dicted position.   The  single  exception  noted  was  for  tra- 
jectory 1.   In  this  specific  case  the  filter  was  surprisingly 
insensitive  to  variations  in  the  elements  of  the  Q(K)  matrix, 

and  its  predicted  position  response  was  fairly  flat  over  a 

-11      -  7 
four-decade  range  (10    -  10   ) .   The  elements  reported 

for  Filter  1A  are  therefore  somewhat  arbitrary,  lying  in 

the  center  of  the  flat  region.   Values  for  the  pair  q^j   =  Q33 

above  and  below  this  region  produced  inferior  results.   A 

refined  search  grid  in  the  flat  region  failed  to  produce 

any  significant  change  in  the  results  reported  for 

Filter  1A. 
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Table  VII  gives  the  simulation  results  for  each  optimal 
filter  configuration  and  for  each  optimal  filter  configuration 
against  the  two  remaining  trajectories.   Table  VIII  compares 
the  hit/scare  percentages  from  the  simulations  of  the  models 
using  the  two  different  state  vectors.   The  data  presented 
for  the  spherical  coordinate  state  vector  represent  the 
optimum  performance  for  the  trajectories  considered.   The 
results  reported  for  the  Cartesian  Coordinate  state  vector 
model  represent  the  best  of  the  results  reported  by  Parr 
[Ref .  4] .   Tables  IX  through  XI   show  the  tracking  per- 
formance for  each  optimal  filter  configuration  against  its 
specific  trajectory.   The  results  shown  are  the  time- 
averaged  means  of  estimation  error  and  the  RMS  value  of  the 
deviation  of  the  means  of  estimation  error  from  their 
time  average. 


42 


to 

<D 

00 

LO 

CM 

U 

<* 

o 

LO 

>H 

rt 

• 

• 

• 

ei 

u 

cr> 

o 

CO 

o 

V) 

T-i 

i-l 

E- 

o\° 

CJ 

W 

»"3 

•p 

LO 

to 

s 

cr> 

o 

CTi 

H 

o\<=> 

oo 

o 

CM 

co 

Q 

J 

W 

w 

CM 

Q 

CM 

0) 

!>. 

CX> 

r-- 

i—i 

O 

*-i 

CM 

OO 

CM 

S 

S 

>H 

cd 

• 

• 

• 

i— i 

Pi 

o 

CO 

«* 

CO 

H 

2 

O 

c/i 

I— 1 

PL, 

o 

E- 

CM> 

o 

i— i 

H 

CJ 
W 

m 

<t* 

•"3 

■P 

o 

to 

VO 

K 

K 

«£ 

•H 

CM 

■*t 

i— 1 

H 

W 

c2 

X 

• 

• 

• 

1— I 

J 

E- 

©\o 

r-^ 

vO 

t-~ 

>— I 

Pi 

w 

CM 

> 

O 

P-. 

u 

(D 

< 

rH 

co 

x> 

H 

s 

0$ 

►J 

o 

H 

CO 

Q 

z 

O 

t-H 

<D 

i>- 

r-~ 

LO 

o 

W 

^ 

CM 

LO 

LO 

(— ( 

H 

>- 

rt 

• 

• 

• 

H 

< 

pi 

u 

t-- 

LO 

r-- 

< 

•-J 

o 

C/) 

t-H 

i-q 

W 

H 

o\o 

p 

Pi 

U 

S 

Pi 

W 

h- ( 

o 

^ 

■p 

rH 

CT> 

en 

CO 

u 

<tj 

•H 

<5t 

CO 

o 

c2 

X 

• 

• 

• 

H 

o\=> 

CM 

to 

Pi 
w 

E- 

i-c 

Pi 


< 

CM 


< 

to 


43 


CO 

,-J 

w 

p 

^ 

o 

< 

HH 

CO 

2: 

w 

o 

H 

»— i 

Pi 

E- 

< 
cj 

P 

►J 

Z 

w 

< 

cj 
cj 

k4 

< 

< 

CJ 

S 

i— ( 

o 

q< 

P 

w 

:z 

33 

< 

P< 

Pi 

(— 1 

CO 

h-i 

n 

1— 1 

W 

m 

> 

X 

H 

H 

< 

<D 

j 

i— 1 

Pi 

w 

,0 

O 

Pi 

a 

Ph 

Pi 

H 

o 

CO 

u 

H 

hJ 

Pi 

p 

o 

CO 

H 

w 

u 

Pi 

> 

Ph 

o 

m 

^ 

< 

o 

H 

CO 

CO 

t-( 

Pi 

W 

< 

H 

PjC 

< 

s 

^ 

o 

(— i 

cj 

Q 

Pi 
o 
o 

CJ 

a> 

J-< 

CO 

W 

rt 

w 

cj 

u 

1—) 

2: 

to 

Pi 
o 

S 

o\« 

H 

i 

cj 

cj 

4-> 

W 

2: 

•H 

►"3 

o 

X 

< 

►J 

o\« 

H 

o 

CD 

w 

f-< 

^ 

cj 

a 

h- 1 

"Z 

u 

Pi 

CO 

w 

cc, 

o\o 

> 

i 

P 

E- 

m 

Pi 

+-> 

2 

o 

•H 

< 

EC 

^ 

s 

CO 

o\° 

0 

f-l 

CO 

w 

rt 

H 

cj 

u 

t— < 

2 

t/> 

Pi 

<tj 

o\° 

o 

Pi 

E- 

i 

CJ 

CJ 

4-> 

w 

z 

•H 

>-D 

o 

^ 

< 

J 

o\o 

Pi 

E-h 

>H 

E-h 

1— f 

O 

O 

-J 

CD 

tq 

w 

fH 

> 

CJ 

Oj 

i 

z 

CJ 

H 

<£ 

CO 

2: 

pc, 

o\«> 

< 

i 

H 

C-H 

CO 

Pi 

•M 

z 

o 

•H 

o 

EC 

^ 

CJ 

CO 

o\<= 

Pi 
o 

H 

u 
w 
> 

w 

H 

< 
E-h 
CO 


P 
W 
E-h 
CO 
W 
E-h 

E-h 
O 

2: 


CXi 
CO 


I— I 


CV1 


CO 

t-l 


CTl 


cm 


<N1 


P 

W 

H 
co 

W 
E-h 

E- 
O 


o 
o 


en 


o 

CM 


o 


o 


vO 


O 
1^ 


o 
•3- 


K) 
f^ 


W 

w 

-J 

H 

:z 

E-h 

< 

< 

< 

< 

cj 

2 

i— i 

^ 

(— 1 

H- 1 

CO 

1— 1 

Pi 

P 

w 

P 

w 

Pi 

H 

Pi 

EC 

o 

ffl 

o 

p. 

o 

< 

o 

CO 

CJ 

CJ 

CJ 

44 


Table  IX 

POSITION  ESTIMATION  PERFORMANCE  OF  THE  OPTIMIZED 

CORRELATED  RANDOM  ACCELERATION  MODEL  FILTERS 

FILTER 

PARAMETER       1A            2A  3A 

e   (yds)      5.87xl0"2     0.111  1.06xl()"2 

o-      (yds)      0.125         0.289  0.135 
er 

e   (deg)      -5.45xl0~5     6.59xl0~5  2.55xl0"4 

a-  (deg)       7.55xl0"4     2.37xl0~3  8.23xl0_4 

e"e  (deg)       1.48xl0"4     2.41xl0-4  1.94xl0"4 

a-  (deg)       7.43xl0"4     1.40xl0"3  7.34xl0'4 
eG 
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PARAMETER 

e^   Cy/s) 

<*£-.  Cy/s) 

r 

e-     (d/s) 

Table  X 
RATE  ESTIMATION  PERFORMANCE  OF  THE  OPTIMIZED 
CORRELATED  RANDOM  ACCELERATION  MODEL  FILTERS 


FILTER 
1A  2A  3A 

3-?7  3.83 

29.71  29.70 


3.78 
29.70 


2.20xl0'4  1.64xl0~2  5.06xl0"3 

1.76X10"3  0.109  1.39xl0"2 

eg  (d/s)      -1.58xl0~4  -l.SOxlO"4  1.63xl0-3 

^.Cd/s)       3.31X10"3  8.14xl0"3  2.67xl0-3 
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Table  XI 

ACCELERATION  ESTIMATION  PERFORMANCE  OF  THE  OPTIMIZED 

CORRELATED  RANDOM  ACCELERATION  MODEL  FILTERS 

FILTER 


PARAMETER 
e^  (y/s 


a-.,  (y/s 
er 


e-    (d/s 


o-..(d/s 
eg  (d/s 


o-   (d/s 

ee 


1A 


0.546 


3.32 


2.68x10 


-9.70x10 


1.52x10 


2A 


-0.549 


3.39 


4.35xl0~4    1.82xl0"2 


0.31 


0.49 


0.10 


3A 

-0.526 
3.17 

-1.12x10 
1.39x10 
5.39x10 
2.67x10 


-4 
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IV.   ADAPTIVE  FILTERING 
A.   THE  ADAPTIVE  BANDWIDTH  FILTER 

A  single  fixed-parameter  filter  configuration  cannot 
realistically  be  expected  to  perform  well  against  all 
trajectory  profiles  likely  to  be  encountered.   One  way  of 
approaching  this  problem  is  to  use  a  filter  that  adapts 
its  characteristics  to  a  changing  tracking  environment, 
e.g.,  a  target  maneuver  is  detected.   The  design  problem 
now  becomes  two-fold;  under  what  circumstances,  or  when, 
should  the  filter's  characteristics  be  changed  and  by  how 
much  should  the  filter  be  changed?  Price  and  Brown  [Ref. 
11]  discuss  some  interesting  approaches  to  the  design  of 
adaptive  filters;  unfortunately  their  complexity  is  beyond 
the  defined  scope  of  this  thesis. 

The  single  method  to  be  considered  here  will  be  a 
simple  s„itch-on-range  scheme.   The  motivation  for  such  a 
scheme  lies  in  the  simplicity  of  its  implementation. 
Additionally,  Parr  reported  [Ref.  4]  that  of  the  methods 
he  investigated,  the  switch-on-range  method  provided  the 
best  overall  performance  in  his  simulations.   When  filter 
tracking  is  initiated,  a  low  Q(K)  gain  schedule  (corre- 
sponding to  a  narrow  bandwidth  filter)  is  used  to  provide 
the  optimum  performance  against  long-range  targets.   When 
the  target's  measured  range  reaches  a  predetermined 
switching  point,  the  filter  switches  to  a  high  Q(K)  gain 
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schedule  (corresponding  to  a  wideband  filter)  for  the 
short-range  tracking  situation. 

Both  filter  models  developed  in  Chapter  II  are  simulated 
in  a  switch-on-range  adaptive  mode  of  operation.   Again, 
to  duplicate  Parr's  testing  environment,  a  switching  range 
of  3333.33  yards  is  used,  and  the  filters  are  tested 
against  trajectories  9  -  14.   The  profiles  for  trajectories 
9-14  are  identical  to  Parr's  trajectories  1-6.   The 
results  reported  for  the  adaptive  simulations  are  the  pro- 
ducts of  80-member  ensemble  Monte  Carlo  simulations.   For 
the  simulations  of  the  constant-velocity  model,  the  low 
Q(K)  gain  schedule  used  was  that  from  Filter  3.   The  data 
in  Table  II  show  that  Filter  3  gave  the  best  overall  per- 
formance against  the  long-range  trajectories  treated  in 
Chapter  III.   Similar  considerations  yield  Filter  4  as 
the  source  for  the  high  Q(K)  gain  schedule  for  the  short- 
range  tracking  environment.   When  filter  tracking  is 
initiated,  the  Kalman  gains  used  are  those  for  Filter  3. 
When  the  target's  measured  range  becomes  less  than  or  equal 
to  the  switching  range  of  3333.33  yards,  the  Kalman  gain 
schedule  is  reset  to  the  steady-state  gain  values  obtained 
from  the  gain  schedule  for  Filter  4.   Similar  consideration 
of  the  data  presented  in  Table  VII  gives  Filter  2A  as  the 
source  of  the  high  Q(K)  gain  schedule  while  Filter  3A 
provides  the  low  Q(K)  gain  schedule.   As  in  the  above 
simulation,  the  low  Q(K)  gain  schedule  is  used  until  the 
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switching  range  (3333.33  yards)  is  reached.   At  that  point, 
the  Kalraan  gain  schedule  is  reset  to  the  steady-state 
gain  values  from  gain  schedule  for  Filter  2A.   Initializa- 
tion of  the  filters  was  the  same  for  all  simulations  con- 
ducted and  is  described  by  Equations  (2.5),  (2.8),  and 
(2.9).   The  statistics  assumed  in  all  simulations  for  the 
measurement  noise  process  v(K)  are  given  by  Equations  (1.6) 
and  (3.1)  . 

B.  SIMULATION  RESULTS  FOR  THE  CONSTANT-VELOCITY  MODEL 
SWITCH- ON- RANGE  ADAPTIVE  FILTER 

The  constant-velocity  model  filter  was  tested  in  a 
switch-on-range  adaptive  mode  of  operation  against  trajec- 
tories 9  -  14.   Table  XII  shows  the  filter's  position  pre- 
diction accuracy  results.   The  data  tabulated  in  Table  XIII 
compare  the  results  for  the  adaptive  simulations  of  the 
models  using  the  two  different  state  vectors.   Tables 
XIV  and  XV  show  the  tracking  performance  of  the  filter  in 
the  adaptive  mode.   As  in  Chapter  III  the  tabulated  data 
represent  the  time-averaged  means  of  estimation  error  and 
the  RMS  value  of  the  deviation  of  the  means  of  estimation 
error  from  their  time  average. 

C.  SIMULATION  RESULTS  FOR  THE  CORRELATED  RANDOM 
ACCELERATION  MODEL  SWITCH-ON-RANGE  ADAPTIVE  FILTER 

The  correlated  random  acceleration  model  was  also 

tested  in  an  adaptive  mode  of  operation  against  trajectories 

9-14.   Table  XVI  shows  this  filter's  position  prediction 

performance  results.   Tables  XVII  -  XIX  show  the  filter's 
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Table  XII 
SWITCH -ON -RANGE  ADAPTIVE  FILTER  SIMULATION 
FOR  THE  CONSTANT-VELOCITY  MODEL 


TRAJECTORY 

9 

10 

11 

12 

13 

14 

SHIT 


90.58 
60.57 
71.58 
63.95 
67.37 
74.67 


I SCARE 

5.55 

12.98 

11.39 

15.45 

9.98 

9.65 
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Table  XIII 
COMPARISON  OF  RESULTS  FOR  THE  SPHERICAL  AND 
CARTESIAN  COORDINATE  STATE  VECTOR  CONSTANT - 
VELOCITY  MODEL  SWITCH- ON- RANGE  ADAPTIVE 

FILTERS 

STATE  VECTOR 


TRAJECTORY 


SPHERICAL 


iHIT    ISCARE 


CARTESIAN 


iHIT    %  SCARE 


9 
10 
11 
12 
13 
14 


90.58 

5.55 

81.80 

11.60 

60.57 

12.98 

47.20 

10.60 

71.58 

11.39 

60.10 

14.20 

63.95 

15.45 

52.30 

13.40 

67.37 

9.98 

60.20 

15.00 

74.67 

9.67 

54.70 

13.10 
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tracking  performance  in  the  adaptive  mode.   As  in  the  pre- 
ceding section,  the  tabulated  data  represent  the  time- 
averaged  means  of  estimation  error  and  the  RMS  value  of  the 
deviation  of  the  means  of  estimation  error  from  their  time 
average . 
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Table  XVI 
SWITCH- ON- RANGE  ADAPTIVE  FILTER  SIMULATION 
FOR  THE  CORRELATED  RANDOM  ACCELERATION  MODEL 

TRAJECTORY  Shit  Iscare 

9  62.27  10.85 

10  46.72  15.20 

11  60.46  10.76 

12  53.34  13.06 

13  56.67  9.36 

14  56.08  11.23 
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V.   SUMMARY  AND  CONCLUSIONS 

A.   SUMMARY  OF  TEST  RESULTS 

1.   Summary  of  Simulation  Results  for  the  Two  Models 
The  results  presented  in  Chapter  III  show  that 
against  non-maneuvering  target  profiles  the  optimized 
constant-velocity  models  provided  superior  prediction 
performance.   When  the  bandwidth  of  the  filter  is  increased, 
adequate  position  prediction  performance  against  maneuvering 
targets  is  obtained  but  at  the  expense  of  performance 
degradation  against  the  non-maneuvering  profiles.   The 
optimized  correlated  random  acceleration  model  showed 
only  a  slight  improvement  against  the  maneuvering  profiles 
considered,  but  produced  results  inferior  to  those  obtained 
from  the  constant-velocity  model  against  the  non-maneuver- 
ing target  profiles.   The  performance  of  the  constant- 
velocity  model  against  long-range  maneuvering  profiles 
was  significantly  lower  than  for  the  short-range  maneuvering 
profiles.   This  is  to  be  expected  since  the  long-range 
profiles  require  that  the  predictions  be  made  over  a  longer 
period  of  time  during  which  the  target  may  maneuver 
unexpectedly  and  significantly.   The  tabulated  data 
presenting  the  tracking  performance  of  the  filters 
developed  in  Chapter  II,  show  that  the  position  estima- 
tion performance  for  both  filters  is  approximately  the 
same.   The  correlated  random  acceleration  model  gave  only 
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slightly  better  rate  estimation  performance  than  did  the 

constant-velocity  model.   Against  the  target  profiles 

considered  the  filter  models  derived  from  a  spherical 

coordinate  state  vector  provided  significantly  better 

overall  predicted  position  results  than  the  filter  models 

derived  from  a  Cartesian  Coordinate  state  vector. 

2 .   Summary  of  Simulation  Results  for  the  Adaptive 
Filters 

The  results  reported  in  Chapter  IV  for  the 
simulation  of  the  switch-on-range  adaptive  filters  indi- 
cated that  the  constant-velocity  model  derived  from  a 
spherical  coordinate  state  vector  provided  significantly 
better  predicted  position  performance  than  the  constant- 
velocity  model  filter  derived  from  a  Cartesian  Coordinate 
state  vector  against  all  target  profiles  considered. 
Further,  the  results  indicate  that  of  the  models  derived 
from  a  spherical  coordinate  state  vector,  the  constant- 
velocity  model  gave  the  best  overall  predicted  position 
performance. 

B.   CONCLUSIONS 

A  filter  derived  from  a  spherical  coordinate  state 
vector  with  assumed  linear  system  dynamics  produces  pre- 
dicted position  results  potentially  superior  to  those 
obtained  from  a  filter  derived  from  a  Cartesian  Coordinate 
state  vector.   The  spherical  coordinate  state  vector 
constant-velocity  model  provides  better  overall  predicted 
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position  results  than  the  correlated  random  acceleration 
model  derived  from  a  spherical  coordinate  state  vector. 
The  reason  for  this  is  perhaps  due  to  the  manner  in  which 
the  predicted  position,  x(K+N/K)  is  generated.   For  both 

.A. 

models  considered,  the  predicted  position  x(K+N/K)  is  formed 
by  raising  the  state  transition  matrix  4>(K,K-1),  to  be 
the  N — -  power  as  shown  by  Equation  (C.2).   In  the  case  of 
the  constant-velocity  model  such  a  method  uses  the  current 
position  and  velocity  estimates  to  generate  the  predicted 
position  x(K+N/K).   The  correlated  random  acceleration 
model  uses  the  current  position,  velocity  and  acceleration 
estimates  to  generate  the  predicted  position.   Thus,  any 
error  in  the  estimate  of  the  target's  acceleration  state 
causes  significant  errors  in  the  predicted  position  when 
propagated  over  the  prediction  interval.   These  accelera- 
tion induced  prediction  errors  are  not  present  when  only 
the  current  estimates  of  position  and  velocity  are  used 
as  in  the  case  of  the  constant-velocity  model. 

A  filter  model  derived  from  a  spherical  coordinate  state 
vector  can  produce  results  potentially  superior  to  those 
obtained  from  the  corresponding  Cartesian  Coordinate  state 
vector  model  in  an  adaptive  filtering  environment. 

C.   SUGGESTIONS  FOR  FURTHER  INVESTIGATION 

Since  the  constant-velocity  model  produced  the  best 
overall  predicted  position  results,  a  complete  analysis 
of  the  sensitivity  of  this  model's  predicted  position 
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performance  to  the  statistics  assumed  for  the  measurement 
noise  process  v(K)  and  the  sampling  period  T  should  be 
carried  out.   Additional  sensitivity  studies  should 
also  focus  on  allowing  the  diagonal  elements  of  Q(K)  to 
vary  independently,  thus  removing  the  constraint  that 

q22  =  q33' 

For  the  correlated  random  acceleration  model,  one  area 

that  warrants  future  investigation  is  the  use  of  only  the 
position  and  velocity  estimates  of  the  filter  in  generating 
the  predicted  position,  and  comparing  the  results  with  those 
obtained  from  the  constant-velocity  model  filter.   This 
would  reduce  the  errors  in  the  predicted  position  due  to 
errors  in  the  current  estimate  of  the  target's  accelera- 
tion state. 

Both  models  should  be  tested  against  a  greater  variety 
of  target  profiles,  especially  targets  that  present  crossing 
profiles  in  which  large  bearing  and  elevation  rates  are 
generated,  and  targets  that  execute  maneuvers  while  in 
level  flight. 

Finally,  more  sophisticated  adaptive  filtering  schemes 
should  be  developed  and  investigated  for  the  constant- 
velocity  model.   A  prime  candidate  for  such  an  investigation 
is  the  residual  testing  method  reported  by  Price  and 
Brown  [Ref .  11] . 
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APPENDIX  A 
MONTE  CARLO  SIMULATION  PROGRAM 

A.l   PROGRAM  LISTING  AND  DESCRIPTION 

The  Monte  Carlo  Simulation  Program  used  in  this  thesis 
was  written  by  Dr .  D.  E.  Kirk  of  the  Naval  Postgraduate 
School.   It  is  extremely  flexible  and  well  documented  to 
assist  the  potential  user  in  data  preparation  and  input, 
as  well  as  the  various  options  available.   It  should  be 
pointed  out  that  SUBROUTINE  KILL  and  all  statements  in  the 
main  program  that  deal  with  miss  distance  data  were  written 
specifically  for  this  thesis.   They  may  be  deleted  with  no 
loss  of  program  flexibility  or  functional  operation. 

Computer  system  requirements  vary  with  the  simulation 
and  the  user's  desires.   All  the  computer  work  done  for 
this  thesis  was  performed  on  the  IBM  SYSTEM  360/67  at  the 
Naval  Postgraduate  School.   The  program  as  currently 
dimensioned  requires  520K  bytes  of  storage.   When  redimen- 
sioned  for  the  60  data  point  trajectories,  320K  bytes  of 
storage  are  required.   A  100-member  ensemble  simulation  of 
a  linear  sixth-order  filter  with  no  on-line  computation 
required,  the  trajectory  read  in,  miss  distance  information 
desired,  and  all  printing  and  plotting  options  invoked 
requires  approximately  1.5  minutes  of  CPU  time  for  execution. 
An  additional  30  seconds  (approximately)  are  required  for 
the  compile/link  process.   For  the  work  done  for  this  thesis 
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the  program  was  pre-compiled  (FORTRAN  IV  Level  G)  and 
stored  and  executed  from  a  data  cell.   The  link  process 
in  this  case  required  approximately  two  seconds.   If  the 
same  linear,  sixth-order  filter  is  again  simulated  but  the 
printed  and  plotted  output  is  suppressed  (retaining  only 
the  miss  distance  data)  59  seconds  of  CPU  time  is  required 
to  execute  the  simulation.   In  the  case  of  the  linear 
ninth-order  filters  that  were  simulated,  the  storage 
requirements  remain  unchanged,  as  do  the  time  require- 
ments for  the  compile/link  process.   A  100 -member  ensemble 
simulation  of  the  linear  ninth-order  filter  with  no  on- 
line computation  required,  miss  distance  information 
desired,  the  trajectory  read  in  and  all  printing  and  plotting 
options  invoked  requires  approximately  3.0  minutes  of  CPU 
time  for  execution.   If  only  the  miss  distance  information 
is  desired  and  all  other  output  is  suppressed,  approxi- 
mately 2.35  minutes  of  CPU  time  is  required  to  execute 
the  simulation. 
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A. 2   LISTING  OF  KEY  VARIABLES  USED  IN  THE  MONTE  CARLO 
PROGRAM 


AMISS 

COVW(m,m) 

EI(n,n) 
ERR(N,nsam) 

G(n,m) 

GAMMA (n,m) 

GAMMAS (n,m) 

GK(n,m) 

GKS (n,m,nsam) 

H(m,n) 

HIT 

HS(m,n) 

ITER 

M 

MISSD(nens ,nsam) 

N 

NENS 
NSAM 
PHI(n,n) 


average  number  of  misses  over  the  entire 
track  (scalar) . 

covariance  of  the  system  random  forcing 
input  matrix  (double  precision) . 

identity  matrix  (double  precision) . 

sample  mean  of  estimation  error  vector 
(single  precision) . 

sample  gain  matrix  (double  precision) . 

system  random  forcing  input  distribution 
matrix  (double  precision) . 

system  random  forcing  input  distribution 
matrix  (single  precision) . 

sample  gain  matrix  (single  precision) . 

filter  gain  schedule  (single  precision) . 

system  observation  matrix  (double 
precision) . 

average  number  of  hits  over  the  entire 
track  (scalar) . 

system  observation  matrix  (single 
precision) . 

ensemble  member  index. 

number  of  measurements . 

sample  miss  distance  matrix  (double 
precision) . 

system  order. 

number  of  members  in  the  ensemble. 

number  of  time  points  to  be  processed. 

state  transition  matrix  (double 
precision) . 
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PHIS(n,n) 
PKK(n,n) 
PKKMl(n,n) 
PKKM10(n,n) 

PKKS(n,n,nsam) 

Q(m,m) 

R(m,m) 

SCARE 

SIGV(m) 

SIGW(m) 

SIGXZ(n) 

V(m) 

VAR(n,n,nsam) 

W(m) 

X(n) 
XHATZ (n) 

XHKK(n) 

XHKKM1 (n) 


state  transition  matrix  (single 
precision) . 

theoretical  covariance  of  estimation 
error  matrix  P(K/K)  (double  precision). 

one-step  prediction  error  covariance 
matrix  P(K+1/K)  (double  precision). 

initial  value  of  the  covariance  of 
estimation  error  matrix  P(0/-1)  (double 
precision) . 

sample  covariance  of  estimation  error 
matrix  P(K/K)  (single  precision). 

system  random  forcing  input  covariance 
matrix  (double  precision) . 

covariance  of  measurement  noise  matrix 
(double  precision) . 

average  number  of  scares  over  the  entire 
track. 

standard  deviation  of  measurement  noise 
vector  (single  precision). 

standard  deviation  of  system  random 
forcing  input  vector  (single  precision) . 

standard  deviation  of  system  initial 
condition  vector  (single  precision) . 

measurement  noise  vector  (single 
precision) . 

actual  covariance  of  estimation  error 
matrix  (single  precision) . 

system  random  forcing  input  vector 
(single  precision) . 

state  vector  (single  precision) . 

initializing  value  of  the  state  estimate 
vector  x(0/-l)  (single  precision). 

state  estimate  vector  x(K/K)  (single 
precision) . 

one-step  state  prediction  vector  x(K/K-l) 
(single  precision) . 
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XM(n,nsam)         mean  state  trajectory  over  the  ensemble 

(single  precision) . 

XS(n,nsam)         reference  trajectory  (single  precision) 

XZMEAN(n)  mean  of  initial  condition  vector 

(single  precision) . 
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A. 3   IBM  SYSTEM  360/67  SOURCE  LIBRARY  SUBROUTINES  CALLED 
BY  THE  MONTE  CARLO  SIMULATION  PROGRAM 


GAUSS3 
OVFLOW 

PLOTP 

SNORM 


matrix  inversion  by  the  method  of  GAUSS 

utility  program  required  by  the  random 
number  generator. 

generates  two-dimensional  line  printer 
plots . 

random  number  generator  that  gives 
normally  distributed,  shuffled  pseudo- 
random numbers. 
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A. 4   MONTE  CARLO  SIMULATION  PROGRAM  FLOWCHART 


START 


XL 


READ 

PROGRAM 

FLAGS 


READ 
<KK,K-1) 


± 


WRITE 
<KK,K-1) 


^ 


116 


117 


READ 


r(K)Q(K)r,(K) 


V 


WRITE 

rCK)QCK)rT(K 

*S»  *"*>  "N* 


READ 
Q(K) 


\/ 


■> 


^  0,1 


<■ 


J^Z 


WRITE 
Q(K) 


118 


t  o 


READ 
P(0/-1) 


± 


WRITE 
PC0/-1) 


t    0 
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© 


± 


\ 


READ 

STD.  DEV.  OF 


M. 


WRITE 

STD.  DEV  OF 
v(K) 


READ 

STD.  DEV.  OF 

. w(K)   


V 


WRITE 
STD.  DEV.  OF 
w(K) 


<r 


120 


READ 
x(0/-l) 


WRITE 
x(0/-l) 


READ 
MEAN  OF 
x(0) 


±. 


READ 
STD.  DEV. 
x(0) 
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0 


0 


WRITE   MEAN 
AND   STD.    DE\4 
OF   x(0) 


*    0 


CALL  "TRACK" 
GENERATE  TRAC* 
BY  SINGLE  CALI 
TO  TRACK 


=  1 


=  0 


\L 


\ 


READ  IN 
SINGLE  TRACK 
TO  BE  USED 


iL 


WRITE  FIRST 
AND  LAST 
POINTS  OF 
TRACK 


<Zr 
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0 


=    1 


CALL  "QMAT" 
GIVEN  Q(K)  FORM 

r(K)Q(K)rT(K) 


<r 


=    -1 


=  0 


V 


CALL  "GAIN" 
COMPUTE  G(K) 
ONCE  BEFORE 
MONTE  CARLO 


^ 


READ 
G(K) 


=  1 


\/ 


> 
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0 


r 


± 

DO 

ITER=1,NENS 


BEGIN  MONTE  CARLO 
LOOP  HERE 


=    0 


<r 


± 


INITIALIZE 
FILTER 


± 


DO    K=1,NSAM 
-CALL    "   MEAS" 
FROM  NOISY 
UF.ASTTKF.MF.NT 


^ 


t    1 


CALL    "XZEROV 
FORM 

x(0) 


<D 
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0 


€> 


=    0 


V 


CALL    "QON" 
COMPUTE 
Q(K)    ON-LIN1 


^ 


=    0 


t  o 


&. 


CALL    "RON" 
COMPUTE 
R(K)    ON-LINIi 


*Zr 


M. 


CALL  "GAIN" 
COMPUTE 
GAINS 
ON-LINE 

> 


© 
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CALL  "TRACK" 

UPDATE 

TRACK 


*    1 


O 


% 

CALL  "ESTIM1 
UPDATE  STATE; 
ESTIMATE 


^ND_MONTE_CARLO    LOOP 


± 


CALL  "STAT" 
UPDATE  SUMS 
FOR 
STATISTICS 


=  0 


CALL  "KILL" 
COMPUTE 
MISS 
DISTANCE 


=  NASM 


f   NSAM 


V 


> 


<r 
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© 


v 

COMPUTE 
STATISTICS  FOR] 
THE  MONTE 
CARLO 

STMTTT.ATTON 


=  0 


=  1 


e 


e 


=  0 


WRITE  G(K) 
P(K/K)  FOR 
EACH  K 


Ak 


WRITE  SAMPLE 
MEANS  AND 
VAR.  OF 
ESTIMATIOI 
ERROfi^FlDR 
K 


PLOT 
G(K)  vs.  K 
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G> 


PLOT 
P(K/K)  vs. 


PLOT  MEAN 


OF 


x(K)  vs  K 


PLOT  SAMPLE 
MEAN  OF 
ESTIMATION 


ERROR  vs  K 


PLOT  SAMPL 
VAR.  OF 
ESTIMATION 


ERROR  vs  K 
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0 


t    0 


± 


=  0 


COMPUTE  MISS 

DISTANCE 

STATISTICS 


f    0 
K — 


<D 


PLOT  MEAN 
OF  MISS 
I STANCE  vs 


K. 


PLOT  STD.  D#V, 
OF  MISS 
DISTANCE  vs 
-K 


PLOT   MEAN 

OF  MISS 

DISTANCE   vs 

ftftN€fi 


Z 


i 


V 


PLOT  STD  DE^ 
OF  MISS 
DISTANCE  vs 
RANGE 
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& 


t  o 


<■ 


/PLOT   MEAN   0] 
-WMISS   DISTANCE 
7vs    PRED.    INWL 


± 


PLOT  STD.  DEV. 
OF  MISS  DISTANCE 
vs  PRED.  INWL 


COMPUTE  FIRING 
DATA  (PERCENT 
HITS,  SCARES, 
HISSES) 


K 


WRITE 

FIRING 

DATA 


WRITE  MEAN 
AND  STD.  DEVI 
OF  MISS  DIS 
TANCE  FOR 
EACH  K 
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APPENDIX  B 
TEST  TRAJECTORIES 

B.l   DESCRIPTION  OF  TEST  TRAJECTORY  PROFILES 

a.  Trajectory  1;  short-range,  low  speed,  non- 
maneuvering 

Trajectory  1  commences  at  a  range  of  3850  yards 

in  a  30  degree  dive  with  an  initial  velocity  of  -230 

yds/sec  (the  negative  velocity  indicates  motion  toward 

the  origin) .   The  trajectory  is  terminated  at  a  final 

range  of  457  yards.   No  maneuver  is  executed.   Sixty 

data  points  were  taken  at  a  4-Hz  rate. 

b .  Trajectory  2;  short-range,  low  speed,  moderate 
level  maneuver 

Trajectory  2  commences  at  a  range  of  3850  yards 

in  a  30  degree  dive  with  an  initial  velocity  of  -230 

2 
yds/sec.   At  3450  yards  the  target  executes  a  2.0  g  sin 

planar  maneuver  lasting  13.5  sec.   The  maneuver  is  completed 

at  the  final  sample  range  of  465  yards.   Sixty  data  points 

were  taken  at  a  4-Hz  rate. 

c .  Trajectory  3;  long-range,  low  speed,  non- 
maneuvering 

Trajectory  3  commences  at  a  range  of  9225  yards 

in  a  30  degree  dive  with  an  initial  velocity  of  -230 

yds/sec.   No  maneuver  is  executed.   The  trajectory  is 

terminated  at  a  final  range  of  5832  yards.   Sixty  data 

points  were  taken  at  a  4-Hz  rate. 
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d.  Trajectory  4;  long-range,  low  speed,  moderate 
level  maneuver 

Trajectory  4  commences  at  a  range  of  9225  yards 

in  a  30  degree  dive  with  an  initial  velocity  of  -230 

2 
yds/sec.   At  8822  yards,  the  target  executes  a  2.0  g  sin 

planar  maneuver  lasting  13.5  sec.   The  maneuver  is  com- 
pleted at  the  final  sample  range  of  5840  yards.   Sixty 
data  points  were  taken  at  a  4-Hz  rate. 

e.  Trajectory  5;  medium-range,  high  speed,  evasive 
level  maneuver 

Trajectory  5  commences  at  a  range  of  4500  yards 

in  a  20  degree  dive  with  an  initial  velocity  of  -300 

2 
yds/sec.   At  3900  yards  the  target  executes  a  6.0  g  sin 

planar  maneuver  lasting  10  sec.   The  maneuver  is  completed 

at  a  range  of  914  yards  and  the  trajectory  is  terminated 

at  a  final  range  of  90  yards.   Sixty  data  points  were 

taken  at  a  4-Hz  rate. 

f .  Trajectory  6;  long-range,  high  speed,  evasive 
level  maneuver 

Trajectory  6  commences  at  a  range  of  8000  yards 

in  a  20  degree  dive  with  an  initial  velocity  of  -300 

2 
yds/sec.   At  7500  yards  the  target  executes  a  6.0  g  sin 

planar  maneuver  lasting  10  sec.   The  maneuver  is  completed 

at  a  range  of  4340  and  the  trajectory  is  terminated  at  a 

range  of  3590  yards.   Sixty  data  points  were  taken  at  a 

4-Hz  rate. 
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g.   Trajectory  7;  medium-range,  high  speed,  non- 
maneuvering 

Trajectory  7  commences  at  a  range  of  5000  yards 

in  a  20  degree  dive  with  an  initial  velocity  of  -300 

yds/sec.   No  maneuver  is  executed.   The  trajectory  is 

terminated  at  a  final  range  of  575  yards.   Sixty  data 

points  were  taken  at  a  4-Hz  rate. 

h.   Trajectory  8;  long-range,  high  speed,  non- 
maneuvering 

Trajectory  8  commences  at  a  range  of  8000  yards 

in  a  20  degree  dive  with  an  initial  velocity  of  -300 

yds/sec.   No  maneuver  is  executed.   The  trajectory  is 

terminated  at  a  final  range  of  3575  yards.   Sixty  data 

points  were  taken  at  a  4-Hz  rate. 

i.   Trajectory  9;  long-range,  low  speed,  non- 
maneuvering 

Trajectory  9  commences  at  a  range  of  8000  yards 

in  20  degree  dive  with  an  initial  velocity  of  -200  yds/sec. 

No  maneuver  is  executed.   The  trajectory  is  terminated 

at  the  final  sample  range  of  50  yards.   One-hundred-sixty 

data  points  were  taken  at  a  4-Hz  rate. 

j .   Trajectory  10;  long-range,  low  speed,  moderate 
level  maneuver 

Trajectory  10  commences  at  a  range  of  8000  yards 

in  a  20  degree  dive  with  an  initial  velocity  of  -200  yds/sec 

2 
At  4000  yards  the  target  executes  a  2.0  g  sin  planar 

maneuver  lasting  13.5  sec.   The  maneuver  is  completed  at  a 

range  of  1200  yards  and  the  trajectory  is  terminated  at 
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the  final  sample  range  of  50  yards.   One-hundred-sixty 

data  points  were  taken  at  a  4-Hz  rate. 

k.   Trajectory  11;  long-range,  low  speed,  moderate 
level  maneuver" 

Trajectory  11  commences  at  a  range  of  8000  yards 

in  a  20  degree  dive  with  an  initial  velocity  of  -200  yds/sec 

2 
At  2700  yards  the  target  executes  a  2.0  g  sin  planar 

maneuver  lasting  13.5  sec.   The  maneuver  is  completed  at 

the  final  sample  range  of  50  yards.   One-hundred-sixty 

data  points  were  taken  at  a  4-Hz  rate. 

1 .   Trajectory  12;  long-range,  low  speed,  moderate 
level  maneuver 

Trajectory  12  commences  at  a  range  of  8000  yards 

in  a  20  degree  dive  with  an  initial  velocity  of  -200  yds/ 

2 
sec.   At  3333  yards  the  target  executes  a  2.0  g  sin 

planar  maneuver  lasting  13.5  sec.   The  maneuver  is  completed 

at  a  range'  of  500  yards  and  the  trajectory  is  terminated 

at  the  final  sample  range  of  50  yards.   One-hundred-sixty 

data  points  were  taken  at  a  4-Hz  rate. 

m.   Trajectory  13;  long-range,  low  speed,  moderate 
level  maneuver" 

Trajectory  13  commences  at  a  range  of  8000  yards 

in  a  20  degree  dive  with  an  initial  velocity  of  -200  yds/ 

2 
sec.   At  2700  yards  the  target  executes  a  2.5  g  sin 

planar  maneuver  lasting  10  sec.   The  maneuver  is  completed 

at  the  final  sample  range  of  50  yards.   One-hundred-sixty 

data  points  were  taken  at  a  4-Hz  rate. 
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n.   Trajectory  14;  long-range,  low  speed,  moderate 
level  maneuver  " 

Trajectory  14  commences  at  a  range  of  8000  yards 

in  a  20  degree  dive  with  an  initial  velocity  of  -200  yds/ 

2 
sec.   At  4000  yards  the  target  executes  a  2.5  g  sin 

planar  maneuver  lasting  10  sec.   The  maneuver  is  completed 

at  a  range  of  600  yards  and  the  trajectory  is  terminated 

at  the  final  sample  range  of  50  yards.   One-hundred-sixty 

data  points  were  taken  at  a  4-Hz  rate. 

B.2   TRAJECTORY  GENERATING  PROGRAM 

The  trajectories  of  Section  B.l  are  the  products  of  a 

tracking  generating  program  written  by  Parr.   A  complete 

2 
description  of  its  use,  development  of  the  planar  sin 

maneuver,  program  listing  and  several  examples  may  be 

found  in  Parr's  work  [Ref.  4]. 
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APPENDIX  C 
THE  N-STEP  PREDICTION  PROBLEM 

The  ultimate  test  of  any  fire  control  system  is  its 
ability  to  destroy  an  engaged  target.   To  this  end,  the 
determining  factor  is  the  accuracy  with  which  it  predicts 
the  future  position  of  the  target,  since  launcher  laying 
orders  are  derived  from  knowledge  of  the  target's  present 
and  future  states.   This  of  course,  neglects  the  response 
characteristics  of  the  launcher  and  the  weapon  (projectile) 
ballistics,  both  of  which  are  assumed  here  to  have 
negligible  effects  on  the  end  result. 

The  interval  over  which  the  predictions  are  made  is 
derived  from  a  straight  line  approximation  of  the  plot  of 
projectile  Time  of  Flight  (TOF)  vs.  the  target's  present 
slant  range,  for  a  5M/54  CAL.  projectile.   For  the  inter- 
val in  question,  {TOF:  TOFe (0 . 0  ,10 . 0) }  the  straight  line 
approximation  takes  on  the  value 

TOF  =  1.43  x  10"3  (sec/yd)  x  (slant  range)    (C.l) 

The  TOF  interval,  (0.0,10.0)  was  based  on  the  assumption 
that  the  maximum  effective  range  of  the  5"/54  CAL  gun 
against  an  air  target  is  7000  yds,  corresponding  to  a 
TOF  of  approximately  10.0  seconds. 
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In  the  absence  of  any  deterministic  inputs  to  the 
system,  the  target's  future  position  predicted  over  N 
samples  is  given  by 

x(K+N/K)  =  cj>Nx(K/K) 


(C2) 


Raising  the  state  transition  matrix  <J>(K,K+1)  to  a  decimal 
power   presents  some  computational  difficulties,  so  N  is 
formed  by  computing  the  TOF  from  x, (K/K)  (the  current 
estimate  of  the  target's  slant  range)  then  truncating  this 
value  to  give  the  integer  N.   For  large  N  the  difficulty 
of  raising  a  (9  x  9)  matrix  to  the  N —  power  still  exists. 
Note,  however,  that  the  assumption  of  linear  plant  dynamics 
decouples  the  filter  in  each  of  three  coordinate  directions 
Partitioning  and  decoupling  we  have  for  the  system, 


x(K+l)  = 


where 


$11     ?        2 


0     0oo         0 
~      ~  z  u         ~ 


33 


xx(K) 


.  ..  -p 
[r,r,r] 


x2(K)  =  [<fr,i,$] 


x3(K) 


[6,e,e]T 


XjCKf 

1l 

0 

0 

w^K) 

x2(K) 

+ 

0 

12 

0 

w    (K) 

~w      J 

x3(K) 

0 

0 

r3 

w^K)  =  [wx(K)] 
w2(K)  =  [w2(K)] 


w3(K)  =  [w3(K)] 


(C.3) 


(C.4) 
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Further,  note  that  for  each  of  the  models  discussed  in 
this  thesis,  each  of  the  <f>..  i  -  1,  2,  3  are  lower  tri- 
angular matrices.   Following  the  method  given  by  Parr 
[Ref.  4],  the  elements  of  the  <J>..  i  =  1,2,3  are  operated 
on  individually  and  the  resulting  <J>  -  -  i  =  1,2,3  are 

adjoined  to  form  the  system  state  transition  matrix 

N 
<f>  (K+N,K)  .   The  intermediate  result  follows  quickly  and  is 

shown  on  the  following  page. 
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N 
Finally,  adjoining  the  ^.jj's  gives 


<f>N(K+N,K) 


N 


hi 

0 

0 

0 

*22 

0 

0 

0 

*33 

(C.6) 


The  result  shown  in  Equation  (C.5)  is  easily  programmed 
and  involves  simple  manipulations  without  the  need  for 
extensive  matrix  multiplication.   Note  that  the  result 
shown  in  Equation  (C.5)  is  equivalent  to  (for  the  single 
dimension) 


r(K+N)  =  r(K)  +  f(K)[(N-K)T]  +  |  r (K) [ (N-K)T] 2     (C.7) 


For  the  case  of  the  constant-velocity  model,  <J>..(1,3), 
4>  -  -  C  2 , 3)  ,  and  $--(3,j)  j  =  1,2,3  are  set  to  zero.  The 
prediction  is  then  easily  found  by  using  Equation  (C.2). 

The  magnitude  of  the  N-step  prediction  error  is 
defined  to  be  the  magnitude  of  a  position  vector,  say 


ERRNET,  drawn  between  the  estimate  at  time  t  =  (K+N)T, 
x(K+N/K),  and  the  reference  trajectory,  x(K+N) 


|ERRNET|   =   |x(K+N)  -  x(K+N/K)| 


(C8) 


Note  that  no  reference  is  made  to  the  direction  of  the 

error  vector  since  the  previously  defined  kill  region  is  a 

sphere  of  constant  radius  centered  at  the  reference 

position  x(K+N) . 
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The  primary  interest  lies  in  the  following  three 
regions 

1)  | | ERRNET | |  <  8.0  yds  hit 

2)  8.0  <  1 | ERRNET | |  <  16.0  yds.         scare 

3)  | | ERRNET | |  >  16.0  yds  miss 

To  avoid  the  necessity  for  coordinate  conversion,  the 
error  vector  is  determined  by  a  simple  rotation  of  the 
polar  plane  followed  by  applications  of  the  law  of  cosines 
Consider  the  situation  depicted  in  Figure  4,  where  points 
a  and  b  represent  x(K+N)  and  x(K+N/K)  (position  only). 
Then  from  Equation  (C.8) 


I | ERRNET J  |  =  | |ab| I  =  | |ba| I  (C.9) 

If  the  polar  plane  in  Figure  4  is  rotated  through  some 
angle  3  about  an  axis  constructed  by  drawing  a  line  through 
the  origin  connecting  the  points  (090,270)  (keeping  r.. 
and  r?  fixed)  Figure  5  is  obtained.   The  angle  of  rotation 


is  defined  as 


3  =  min(6  ,92)  (CIO) 


Point  now  lies  on  the  rotated  polar  plane,  the  position 
of  point  b  remains  fixed,  and  point  c  represents  the  fixed 
range  projection  of  point  b  onto  the  polar  plane.   Some 
rather  obvious  geometric  relationships  follow  immediately; 
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AG  =  /boc  =  |e1  -  e2| 

AR  =  \ai\    =    |r1  -  r2|  (C.llj 

I  od  |  =  maxfr^r  ) 

The  problem  then  remains  of  finding  | ab | .   Figure  6  shows 
the  projection  onto  the  rotated  polar  plane.   Note  that 
triangle  cod  is  an  isosceles  triangle,  then 


/cd0  .  180  -  Af  (C.12) 


Applying  the  law  of  cosines  |ac|  is  given  by 

|ac|  =  {(|^|)2  +  [(2.0)(sinM7T)(|o1|)]2  - 


27U 


.M_._-»rn^n  ^./^^  \1/2 


(C.13) 


(2.0)  ( I  ad  I  )  (2.0)  (siny^-)  (|od|)  cos/cdo)  } 


be |  is  given  by 


|  be  |  =  (2.0)  (|od|)  (sin— -)  (C14) 


Again  applying  the  law  of  cosines  and  noting  that 
/bed  =  /bco ,  | ab |  is  given  by 

fab |  =  {|bcV  +  |a^|2  -  (2 . 0)  ( |bc | ) ( | ac | )cos (/bed) }1/2   (C.15) 

Thus,  the  magnitude  of  the  N-step  prediction  error  has  been 
determined  without  the  need  for  any  coordinate  conversion. 
The  validity  and  accuracy  of  this  means  of  determining  the 
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error  magnitude  was  verified  by  conversion  of  x(K+N)  and 
and  x(K+N/K)  to  the  Cartesian  Coordinate  system  and  deter- 
mining the  magnitude  of  the  error  vector  by  standard  means 
In  every  case  investigated,  the  results  given  by  the  two 
different  methods  were  in  complete  agreement.   Again  note 
that  no  attempt  is  made  to  determine  the  direction  of  the 


error  vector  |ab|  since  the  regions  of  interest  specify 
concentric  spheres  of  constant  radii  centered  at  ^x(K+N) . 
The  directionality  of  the  error  vector  is  eliminated  in 
SUBROUTINE  KILL  through  the  use  of  the  absolute  values  of 
the  differences  in  bearing  and  elevation,  and  by  using  the 
maximum  value  of  the  actual  and  estimated  ranges. 
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